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The Sudden Approximation (SA) for scattering of atoms from surfaces is generalized to allow for 
double collision events and scattering from time-dependent quantum liquid surfaces. The resulting 
new schemes retain the simplicity of the original SA, while requiring little extra computational ef- 
fort. The results suggest that inert atom (and in particular He) scattering can be used profitably to 
study hitherto unexplored forms of complex surface disorder. 
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I. INTRODUCTION 



Structurally disordered surfaces have been a subject of great interest for some time now. Of special interest are epi- 
taxially grown films, liquid surfaces, and amorphous surfaces. In epitaxial growth for example, metal or semiconductor 
atoms are adsorbed on a surface under thermal conditions, to form two- and three-dimensional structures on top of it. 
The physical and chemical properties are determined by the final form of these structures. These may be of dramatic 
importance, e.g, in the production of electronic devices. One of the most exciting aspects of epitaxial growth kinetics, 
is that it creates disordered structures in the intermediate stages. The disorder manifests itself in the formation of 
various types of clusters or diffusion limited aggregates on top of the surface. These structures may be monolayers 
(usually at high temperatures, when the diffusivity is large, or at coverages significantly below one monolayer), in 
l— — '■ which case the disorder is two-dimensional, or they may be composed of several layers, giving rise to disorder in three 
I ' dimensions. Epitaxially grown structures of this type offer an exceptional opportunity for both experimental and 
, theoretical studies of disorder. No satisfactory and comprehensive theory of the epitaxial growth process is as of yet 
available, much due to the absence of reliable interaction potentials for the system. The situation with respect to 
C*~) \ liquid and amorphous surfaces is similar: very little is known at this point about their structure. Progress at this 
stage thus hinges critically on data available from experiments. An important experimental technique is thermal atom 
scattering, and in particular He scattering |]l]-^). The main advantage offered by He scattering is complete surface 
sensitivity, as thermal He atoms do not penetrate into the bulk, unlike other scattering techniques such as neutron 
or X-ray scattering, or LEED. Another important advantage is that He scattering is highly non-intrusive, due to the 
inertness and low energy of the He atoms. The latter also means that He scattering is really a diffraction experiment at 
the typical meV energy scale at which most experiments are performed, with sensitivity to atomic-scale features. The 
interpretation of He scattering experiments is, however, rather involved due to the complicated interaction between 
i ■ the He atom and the surface. 

A highly successful theoretical method in the study of He scattering is the Sudden Approximation (SA), introduced 
in the context of atom-surface scattering by Gerber, Yinnon and Murrell in 1978 ||. It is to date one of the most 
useful tools in the field, and has been reviewed by Gerber |||J. While the original formulation of the SA applied to 
periodic surfaces, it was extended by Gersten et al. to deal with scattering from disordered surfaces |1G[ |. The SA was 
subsequently used to study a large variety of disordered systems, such as isolated adsorbates on crystalline surfaces 
jy] Q, mixed overlayers of Xe+ Ar Jl5|Jl^| and Xe+Kr on Pt(lll) randomly corrugated hard walls j[7), vacancies 
and CO adsorbates on Pt(lll) ]18| , |19| , percolation lattices of substitutionally disordered Xe+Kr monolayers p0[ , 
amorphous mixed monolayers and liquids [glj , compact islands and diffusion limited aggregates , and fractals 

p5| , p6[ . Common to many of these studies was the achievement of detailed insight into the collision dynamics, the 
prediction and understanding of interesting features in the angular intensity distribution of the scattered atoms. In 
several cases it was shown that the angular intensity distribution exhibits nonspecular maxima of two types: Several of 
the peaks are defect-induced rainbows effects, while others (at angles nearer to the specular) are Fraunhofer diffraction 
interferences. Both types of peaks contain useful, largely complementary, information on defect geometry and on the 
Hc/dcfcct interaction. Several further extensions of the SA include the study of Schinke and Gerber which treated 
phonons and the resulting Debye- Waller attenuation |27|; and Hinch's work |^8| which extended the SA to encompass 
larger parallel momentum transfers. 
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The main purpose of this work is to extend the SA in two useful directions: the inclusion of double collision 
processes, and scattering from time-dependent quantum liquid surfaces. Unfortunately, at the time of writing exper- 
iments are unavailable for comparison with the results obtained here. The developments will therefore be primarily 
methodological, in anticipation of experimental data. It is hoped that the results obtained here will motivate He 
scattering experiments on disordered solid and liquid surfaces. As demonstrated by this work and others before it, 
He scattering can provide a wealth of information on disordered surface structure and dynamics. 

For reasons detailed in Sees. || and EH, the SA in its original formulation cannot include double collision events. 



These, however, are important in many cases, especially when the surface corrugation is large. An analysis of 
"rainbow" peaks arising due to double collisions and performed here for the first time, shows that these are distinct 
from single-collision rainbows in that they have an incidence-energy dependence. This should be clearly observable 
experimentally, e.g., in scattering from CO defects on Pt(lll). Multiple scattering by defects has been considered 
in the past, e.g., by Jonsson, Weare and Levi p9| ] who derived interesting results concerning the dependence of the 
specular intensity on coverage for different models of disorder (applied to Xe and CO); by Kara and Armand p0| , 
who considered elastic scattering from a Cu(llO) surface covered with randomly distributed copper atoms; A similar 
multiple scattering theory was presented by Armand and Salanon fjlf , who showed that for randomly distributed 
defects, the incoherent cross section is equal to the product of a form factor and a Fourier transform of the two 
body correlation function of the defects' positions. Common to these works is a multiple scattering expansion of the 
Lipmann-Schwinger equation for the T-matrix, in powers of the number of He scattering events E9], a perturbation 



potential due to defects 30 1, or the number of defects involved in the scattering |3l). The approach to be presented 
here is somewhat similar to that of | p9| , in that the He atom is assumed to undergo a second collision. However, the 
combination of the SA and the Born approximation used here is new, and leads to expressions of significantly greater 
computational simplicity than that of the general multiple scattering theory. 

The study of liquid surfaces is an exciting and novel topic in its own right []32|-[34||, and to date there are very 
few theoretical methods available for this purpose, in which the surface degrees of freedom are treated classically 
P,R5|. The extension of the SA that is presented here, in Sec.tV, treats the surfaces as a quantum liquid, within a 



time-dependent self-consistent-field (TDSCF) approach. This methodological part of the paper serves to generalize 
the SA to include time-dependent processes and thus renders this useful method significantly more widely applicable. 
Concluding remarks are brought in Secjv|. 

II. BRIEF REVIEW OF THE ORIGINAL SUDDEN APPROXIMATION 

Consider a He atom with mass [i incident upon a surface with wavevector k= (K, k z ). KK. — h(k x , k y ) and KK' are 
respectively the intial and final momentum components parallel to the surface, and hk z is the incident momentum 
normal to the surface. The position of the He atom is r = (R, z), where R = (x,y) is the lateral position. The SA 
is valid when the collisional momentum transfer q = K' — K in the direction parallel to the surface is much smaller 
than the momentum transfer normal to the surface: 2k z S> |q|. This condition is satisfied close to specular scattering 
and energy E — (Tik) 2 / (2m) , and moderate surface corrugations. When it holds, one can approximately consider the 
scattering along z as occurring at fixed R. Then if tp is the He wavefunction, it satisfies a Schrodinger equation where 
the dependence on R is adiabatic: 

~^+V R (z)^ R (z) = EMz). (1) 

Here V R (z) is the He-surface interaction potential and no inelastic channels are included. This means that each 
surface point R gives rise to an elastic real phase shift ?7(R), which can be evaluated in the WKB approximation from 
Eq.(§ as: 



r?(R) = / dz 



2m, 



k/ - -2"U R (2) 

Ti 



1/2 



fe^(R), (2) 



where £(R) is the classical turning point pertaining to the integrand in Eq.(Q). The phase shift in turn yields the 
S-matrix as: <S(R) = exp[2i?/(R)]. The R coordinate is conserved in this picture so the S-matrix is diagonal in the 
coordinate representation: 

(R'|5|R) = e 22 " (R) ( 5(R' - R). 

Experimentally one measures probabilities |5(K — > K')| 2 for K — » K' transitions. To obtain these (R|K) =exp(iK- 
R) / \[A (where A is the area of the surface) can be used, to find: 
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(K'|«S|K) = f dR' dR(K'|R')(R'|5|R)(R|K) = - J dRe'^e 2 ^^ . (3) 
This is the well-known expression for the SA scattering amplitude p|jL0| . 

III. THE ITERATED SUDDEN APPROXIMATION FOR DOUBLE COLLISION EVENTS 

A. Background 

The SA developed in Secjfi] can inherently describe only single collision events between the incident particle and 
the surface. This is because the SA assumes that the momentum transfer in the direction parallel to the surface is 
much smaller than in the direction perpendicular to it, i.e., the surface is assumed "not too corrugated". A more 
formal way to see this is by realizing that the SA essentially treats the R = (x, y) coordinate as adiabatic, which 
implies that each trajectory takes place at constant R. Clearly, no double collisions can occur under such conditions. 
However, double collisions may take place, e.g., when an incident atom is scattered off a defect onto the surface, or 
in the opposite order, as depicted in Fig. |[ The purpose of this section is to develop an "Iterated Sudden" (IS) 
approach, which can be interpreted to include double collision effects as well. 

Starting from the first principles of non-relativistic scattering theory, the exact expression for the T-matrix element 
is @: 

t(k-k') = {k'|V|k+) 1 (4) 

where |k+) = f2 + |k). Here |k) is the incident state (£ = — oo) and f2 + is the M0ller operator, which propagates |k) to 
the state at the time of collision at t = 0. In the Born approximation it is assumed that |k+) « |k) |3(|, which yields 
f(k— >k') = (k'|V|k). The IS, on the other hand, consists in setting: 

|k+) » |k sud ) ; |k sud ) = S sud |k), (5) 

and where S sud is the S-matrix operator in the SA [Eq.([|)]. Hence: 

t(k -> k') m i sud (k -> k') = (k'|y|k sud ) (6) 

Eq.(||) can be argued to describe a double collision process as follows. In the exact expression Eq.(||), |k+) represents 
the actual state at t = ("the moment of collision"), as it evolved from t — — oo. However, |k sud ) is the state at 
t = +ooJ^] in the SA. Thus Eq.(||) implies setting: 

"actual state at collision" — "approximate state much after single collision". 

In other words, at the moment of collision, the wave function is taken to have already undergone a prior collision. 

A qualitative argument can be given to justify the validity of this approximation. The Lipmann-Schwinger equation 
for the exact wave function ip is: 

ij = $ in + G Vi>. (7) 

Now, an equivalent way to view the IS is by setting ip w tp sud in the right-hand side of Eq.(^): 

= [<S> in + G Vr nd ] + G Q VAi>, 

where Atp = ip — tp sud . Thus by assuming tp w ip sud , Gq V Aip is neglected. To estimate the magnitude of this term, 
consider two length-scales: First, it is clear that Go V Atp is small at large distances from the surface, since V becomes 
small. Secondly, Go V Atp is also small close to the surface, due to the nature of the SA: the SA diagonalizes the 
Schrodinger equation (T + V)ip = Etp~, E^j = en Sij by taking the energy matrix E to be a scalar matrix: E — eo I, 
6q = m&XiCu ||]. Hence, the SA essentially assumes that the other (i > 0) diagonal elements €4 are negligible in 
comparison to the elements of the diagonalized V. This condition is well satisfied close to the surface, where V is 
largest. Hence one expects Atp « to good accuracy for both large and small distances from the surface. Interpolating, 
one can reasonably expect Gq V Aip to be negligible at all distances. 



1 The M0ller operator fi+ propagates the in-state \ipi n ) to the actual state at t = 0. The S operator propagates IV'in) to the 
actual state at t — +00 pql. 
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B. Calculation of the T-Matrix Elements 

The calculations in this section will be performed in two ways: the first will lead to a result which can be used in 
numerical simulations, the second will employ some simplifications in order to cast the results in a form amenable to 
analytic analysis. 



1. Wave Function Iterated Sudden 

From the close-coupling expansion it follows [fio| that the Sudden wave function in the continuum case is: 

|<(r)) = J dq|^(r)>, 

where j = 1,2 depending on whether one is in the classically allowed or forbidden region, and: 

(r|V q (r)> =$^) e i(K t n+q)-R 



The last equation is where the Sudden approximation is made: Vf n ( r ) was replaced by <&^(z), i.e, the R coordinate 
is treated as an adiabatic parameter. This is now further evaluated in the WKB approximation: 



■. exp 



/ dz' fc(R, z') 



Here: 



fc(r) = Jk 



In the classically allowed (forbidden) region fit 1 ' — —tt/4 and fc(r) is real (t^ 2 - 1 = 0; fc(r) is pure imaginary). The 
normalization ((ipf n \^Pi n ) = 1) is found after some algebra: 



a 



dv |<F R (z) 



(27T) 



dr 



k(r) 



Now it is possible to evaluate the matrix element within the IS: 

**(k5n - 1W) = (kout|^(r)|<) = 

J dq J dvdv' (k out |r) (r|F(r)|r') (r'|^(r)) = 

J dq J drdr' (2tt)~5 e -^t-r [y( r ') 5( r _ r ')] e ^K in +q)-R' = 

— [ drV(r)e^ AK - R + k ^ z ^ [ dB.' $Uz) fdqe^ (R - R '^ = 
(2tt) 2 J J J 

V2^a J drV(r) e "^K-R + k* ut z) ^(z). 

where AK = K out — K; n . Employing the WKB approximation for the wave- function ^^(z), one has finally for the 
transition amplitude from the disordered surface within the IS: 



(8) 



^(k in 



a e 



2tt / dr 



V(r) exp 



( / dz' jfc(R, z') + AK • R + k z out z ) 



(9) 



Comparing the IS expression (||) to the SA expression (^), note that the computational complexity is essentially 
increased merely by an additional z-integration. Considering the computational simplicity of the SA, the IS of fers an 
attractive improvement. A further discussion of the physical significance of the IS result is given in Appendix A 1. 
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2. Phase-Shift Iterated Sudden 



Equation (^) is suitable for numerical applications. In order to make further analytic progress, I now develop an 
alternative, more approximate expression. This will involve calculating the T-matrix element for the IS again, but 
assuming that o/f-shell (i.e., non energy conserving: |q| 2 7^ |q'| 2 ) S-matrix elements are well approximated by on-shell 
ones. The on-shell elements to be used are just the SA elements from Eq.(|J). That these are indeed on-shell, follows 
simply from the fact that Eq.(|) represents the SA result for elastic scattering. Starting again from Eq.(^) [standard 
intermediate steps are omitted from now on for brevity; the details are as in Eq.(0)]: 



t»'(k ta -» k out ) = (k out |F(r)5 sud |k in ) = 



-i (k ut* r — &f n z ■ 



R -^s q , q y(r) 



Due to the double integration over the unrestricted intermediate momentum values q, q', the S-matrix element is 
off-shell. Approximating Sq>q by the on-shell expression Eq.(||) yields: 



— ^/ dq [J d RV(K,„-q)-R' e 2^(R')^ J dry(r)eJ [Rq + fcf„ Z -k out . r] 



(27T) 



And finally: 



i J (k 



in * ■Kout J 



(27T) 



t)A J 



(10) 



(11) 



The last result is the phase-shift equivalent of E q . (^[) , in the on-shell approximation. This form lends itself more 
easily to analysis, as will be demonstrated in Sec. Ill C . Appendix A 2 provides a further discussion of the physical 
significance of this result. 



C. Analysis 



In this section I will analyze the expression for the scattering amplitude, including double collisions [Eq.(ll)], by 
considering a simple model potential, of the form: 

V = V[z-£(R)]. 

This is a rather general potential, in that the surface shape function £(R) can be chosen arbitrarily More examples 
can be found in Ref. |57]]. The significance of the V[z — £(R)] potential is the following: z — £(R) is a two-dimensional 
surface in real space, and is clearly also an equipotential surface. Thus z — zq + £(R) defines a continuum of 
equipotentials as a function of the height zq above the surface plane. Hence V[z — £(R)] can be interpreted as the 
potential of a corrugated surface, with £(R) representing the corrugation. Now, inserting V into Eq.(|ll|) one obtains: 

t(k in -> k^) oc J dve 1 < 2 "( R > - Ak ' r > V[z - e(R)]. 

Transforming to 9 = z — £(R) yields: 

t(k in ^k out )«| dxdyd9V{6)e l ^)-^K-^(^) + e)] = 

y"dRe ! [ 2 "( R )- AK - R - Afc ^( R )] J d9V(d)e- lAk * 9 . (12) 

In 

Clearly, only /r, can account for double collisions. To proceed, it is convenient to perform a stationary phase (SP) 
approximation |38|], with the purpose of studying the singularities in the spectrum, i.e., the "rainbow" peaks. The SP 
yields: 
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AK = 2V?7(Ro) - Ak z V£(R ), 



(13) 



and: 



|/ R (AK)| 2 « (2^) 2 



J <9R 2 



Afc ; 



5 2 C(Ro) 



<9R 2 



This singularity is a result of the crude SP approximation, and is smoothed out in the full scmiclassical expression 
Eq. (|l2"|) . However, it provides the condition from which the inflexion point Rq is found: 



, a 2 r?(Ro) 
' <9R 2 



Ak : 



a 2 c(R ) 



<9R 2 



0. 



(14) 



Eqs.([l3|),(|l4|) together define the double- collision rainbow (DCR) condition. After solving for R from Eq.(|l4]), Eq.(|T3]) 
determines the position of the DCR peak. It is similar to the SA rainbow condition |22]| , AK = 2V?7(Ro). The 
difference is thus in the second term, with the interesting feature that it depends on the incidence momentum in 
the normal direction. A direct prediction of the present analysis is therefore that a DCR peak is distinguished 
from a single- collision rainbow peak by a positional dependence on incidence energy. This effect should have a clear 
experimental signature. Further, the actual position of the peak depends on the corrugation of the surface shape 
function at the inflexion point. This position thus provides valuable information on a surface structural parameter. A 
more complete analysis, including numerical calculations of the double-collision rainbow intensities, is deferred to a 
future publication. 



IV. TDSCF SUDDEN APPROXIMATION FOR ATOM-SURFACE SCATTERING 

In this section I will present a generalized, full derivation of the SA. The derivation is based on the original SA 
paper and its extension to disordered surfaces ju^, but instead of considering the surface atoms as static, they 
are assumed to be moving under the influence of a time-dependent self consistent field (TDSCF). The result will be a 
TDSCF-SA, which is presented here for the first time. Applications may include scattering from quantum vibrating 
solid and liquid surfaces. A similar approach was pioneered by Gerber p5[ , but with the surfaces degrees of freedom 
treated classically. A variation on this mixed quantum-classical approach can be found in Ref. Isl . 



A. Preliminaries 



Consider a set of N surface particles moving in a TDSCF, u„ = (x n ,y n ) being the position of the n th particle, 
assumed constrained to the surface z = 0. Let f Kn (u n ,t) = (u rl \K n ) be the spatial wavefunction for a surface particle 
with wavevector K n = (k™, k™), solving the Schrodinger equation 

H n f Kn (u n ,t) = ih — , 

ot 

with the SCF Hamiltonian H n — v n (u n ,t)+T n , where v n is the SCF potential for the n th particle and T„ its kinetic 
operator. A TDSCF separation based on the atoms as the separable degrees of freedom may not be the best, as the 
residual (their interaction) may become large due to collisions. Nevertheless, for a liquid (the main application I have 
in mind) no better separable coordinates are known. Moreover, it is possible to switch to lattice modes (phonons) at 
a later point in the formalism (see below). Now, from the completeness relation for |u n ) and \K n ), 



du n |u„)(u„| 

i t follows that fn n (u n ,t) is orthogonal in both variables: 

J du n f* n (u n ,t)f K : n (u n ,t) = (k„|k4) = 5(K n - k'J (15) 

J dK n /*„(u„,f) f Kn (u' n ,t) = (U n \u' n ) = S(u n - W n ). (16) 

In addition the SCF separation of the surface particles implies: 



(u n |iv) = (k„|k„') = <W. 
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B. Matrix-Diagonalization Sudden (MDS) 



Turning attention to the He particle, described by the wavefunction ip, the time-dependent Schrodinger equation 
reads: 



[T + V(R, z; {u n } N ; t)ty(r; {u n } N ; t) = ih 



dip(r;{u n } N ;t) 

at 



(17) 



where T = (?i 2 /2/i)V 2 the kinetic energy operator of the He particle, V its potential energy which depends on {u„}jv, 
the positions of the TV surface particles. The surface is assumed to be almost frozen during the collision, so that 
the dependence of V on t is very slow. This time-dependence arises from the slow motion of the surface particles, 
whose dynamics I wish to include presently, in order to develop a description which accounts for the details of the 
energy transfer. It should be noted that it is inconceivable at this point to imagine an experiment which is capable of 
resolving the details of the energy transfer at the level of surface atoms in a liquid, which is thermalized very quickly. 
However, such experiments may become possible in the future (indeed, efforts are underway in the context of atomic 
Ar clusters |}9|) so that a formalism capable of describing them is desirable. 

V and the SCF potentials are not independent; assuming that the influence of the He particle on the surface is 
small and localized, it should be sufficient to evaluate V at r = 0, thus avoiding the need to find ip in order to calculate 
the SCF potentials: 



v(u n ,t) = J Jj£iUi/i t (ui,t)V(0;Ui;t)/ ( « 4 (u i ,t). 



(18) 



It is of course possible to average over ip as well; then Eq.([18|) should be regarded as the first iteration in a TDSCF 
scheme which includes the He particle at later iterations. I now assume that the SCF potentials have been obtained 
and the problem is the solution of the Schrodinger equation (^) . For this purpose it is convenient to expand ip in a 
generalized time-dependent close-coupling form: 



N 

Mr;{Un} N ;t) = J2 Jdejdqj dn n e^+ K ) R / Kn (u„, i)$ eqK „ (z) e -' 



-i(E+e)t/h 



(19) 



Here hh = 7i(K, k z ) is the momentum of the incident He particle, Ti(\ — Ti{q x , q y ) the momentum transfer parallel to 
the surface, E= (frk) 2 /(2/.i) is the He incidence energy, and e is the energy transfer between the He and the surface. 
Since periodicity is not assumed, the q's are unrestricted and continuous. Next the potential is taken to have the 
form:R 



z,u n ,t) 

n=l 

L 

V n (TL, z, u„, t)=J2 V nl (z, t)Q nl (R, u„, t). 



(20) 
(21) 



1=1 



No restrictions are placed at this point on the functions V nl and Q nl or on their number L, so the expansion is general 
(more so than in the other derivations of the SA, Refs. |^,|l^]). It is however assumed, in the spirit of SCF, that the 
He atom has only two-body interactions with the surface particles. The close-coupling equations will follow from 
inserting these expressions for ip and V into the Schrodinger equation ([i"7|). Since the surface particles are decoupled 
by virtue of the SCF approximation, this results in N equations of the form: 



de / dn r . 



a 2 



+ 



d 2 



/I 



+ 



2^ \dz 2 9R 2 j 2m n du 2 t 



£ V nl (z,t)Q nl (R,u n ,t) + (E + e) 



f Kn (u n ,t)e-^ E +^ h J dqe^ +K ) R $ eqK „(z) = 0, 



(22) 



2 n, I are indices, not powers. 
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where the term df Kn (u n ,t)/dt was neglected with respect to (E+e)/fr f Kn (u„, t), in accordance with the assumption 
of slow surface dynamics. 

Multiplying throughout by (/i/irh 2 ) exp(— iq'-R) and integrating over R yields, using J dll exp(iq ■ R) — 2nd(q): 



<i< ( / dK n ^2 - (q + K ) 2 

L 

$> n '(z)Q£ q ,(u„) 

,i=i 



ft d 2 



f Kn {u n ,t)$ 
/ K „(u n) t)$ eq ^ n (z) j e-^+^/B = o. 



Here: 



Q£ q ,(u„) = i- y dRe l ^'-^ R Q„KR,u„). 



(23) 

(24) 
(25) 
(26) 



Next the orthogonality of the surface particles' wavefunctions is employed to eliminate the integrals over n n in Eq.(|23j 
Multiplying throughout by /*, (u„, t) ex.p(i(E + e')t/%) and integrating over u„ and t yields, using Eq.(|l5|): 



^-(q + K) 2 + fc £ 2 



$eq K „(z) + J dK' n F? Knte , < $ e , clK , n (z) = 



(27) 



Here: 



Fit 



rn 



- / dtdu n f*, (u n ,t)-^f Kn (u n ,t)e 

n J n OUn 



i(t'-t)t/h 



is (up to a constant) a component of the kinetic energy matrix of the n th surface particle associated with an e — > e' 
transition of the He particle, and: 



Q£L,e' q '< = / rfidu„/*,(u„,i)Q^ q ,(u„)/ K „(u n ,t)e 



i{e'-e)t/h 



(28) 
(29) 



is a component of the matrix of the Fourier-transformed interaction between the He atom and the n th surface particle. 
Defining the diagonal "energy operator" g 2 with non-zero components 



9, 



,2 = ,2 
eq — ft e 



(q + K) 2 , 



(30) 



Eq.(p7|) can be rewritten as follows: 



J de' dq' dn' n 



5[(e-e')/h]5(q-q')6(K n -K' n 



Q z 2 + fe'q' 



5(q-q')F { 



(31) 



This can be thought of as the (e, q, K n ) component of a linear operator equation for the vector |$"), whose components 
are the wavefunctions $ eqKii (z). Clearly then, the He kinetic energy operator on the first line of Eq.(|l]) is fully 
diagonal, the n th surface particle kinetic energy operator T n on the second line is diagonal in q, and the interaction 
operator Ti n on the third line is not diagonal. In the Born approximation one essentially retains only the fully diagonal 
(kinetic) term fl36|. Ideally one would like to diagonalize the complete operator acting on |<fr n ), which would solve the 
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scattering problem. This is of course unrealistic, and here I will settle for diagonalizing the interaction operator, 
is to this end that one introduces the Sudden approximation, which is usually written as fc| : 



It 



2M Z > 7i|q| 



(32) 



This implies that the incident momentum in the direction perpendicular to the surface (hk z ) is large compared with 
the momentum transfer for all the channels that significantly contribute to scattering. Here it is necessary to add the 
assumption: 



£> e. 



(33) 



These conditions are satisfied for most present-day experimental systems. Using Eqs. (|32|) , (|33|) one can neglect all but 
the specular term in the energy operator, i.e., assume that: 



Se ~ din 1 = k 2 z X, 



(34) 



where X is the identity operator, 
somewhat stronger than Eq.(p3): 



Since g 1 = k\ — (q+2K) • q+2/ie/h 2 , it is seen that the real assumption made is 



fc z 2 »|(q+2K)-q|. 



(35) 



This implies that the incident beam should be nearly perpendicular to the surface, a condition which is in fact 
sometimes not met in practice, where the beam is often incident at 45° and the SA calculations have to be corrected 
accordingly [p3| , p4|j4C| l . The utility of the approximation Eq. (|34]) is that it will leave the energy operator g 2 diagonal 
after diagonalization of the interaction operator in Eq. (|3~ij) . However, whereas this suffices in the stationary surface 
case, here one is facing also the operator J 7 ™, which is diagonal in q but not in K n and e. For this reason it cannot 
be diagonalized simultaneously with H n , since to diagonalize this operator it will be necessary to rotate the q basis 
and hence introduce q-off-diagonal elements into T n . Nevertheless the problem can be solved if one assumes as 
before when df Kn (u n ,t)/dt was neglected, in the spirit of the SA, that the kinetic energy associated with the surface 
particles' motion is very small compared to that of the incident He atom. This is certainly correct for solids, and 
also a reasonable assumption for liquids at common experimental temperatures. Then it is possible to simply neglect 
F" Kn e , K , with respect to k 2 z . With these approximations Eq.(|3l|) becomes, in operator form: 



dz 2 



0. 



(36) 



This represents a set of coupled equations for scattering of the He atom by the n th surface particle into the available 
channels indexed by e, q and K n . To uncouple the set one needs a unitary matrix B n that diagonalizes all the (z- 
independent) Q nl matrices [whose components are defined in Eq.(|29|)]. That a single B n can diagonalize all L Q nl 
matrices is not a priori obvious. I will assume for the moment that this is true and come back to it later [see Eq.(|42"|)]. 
Then: 



with matrix elements 



= B n Q nl (B n ) 



(e'q'<|Q2>Wn) = (eqKnlQjo |eq«n)<S(e - e')S(q - q')5(K n - < 



(37) 



Since B n leaves the energy matrix diagonal in view of Eq.(| 
after application of B n Eq.(Bq) becomes: 



and commutes with purely z-dependent operators, 



O 2 



i=i 



13") = 0, 



(38) 



where 



The system of equations (38) is uncoupled since the operator acting on |$") is now diagonal. Physically, this corre- 
sponds to the fact that there are no transitions between the decoupled channels labeled by e, q and K n in the new 
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system. Thus the scattering matrix Sd that arises from Eq.(|38|) is diagonal and represents elastic scattering in each 
of the decoupled channels. Its elements therefore satisfy: 

(e'q<|«SS|eqK„)(t) = e 2 ^W*)<J(e - e')S(q - q')£("n - <), (39) 

where ??™ qK;n (i) is the elastic scattering phase shift determined for the tc\K n channel from Eq. ( ]38| ) . The time-dependence 
of the S-matrix element may appear strange (after all the S-matrix describes the entire scattering process from 
t = — oo — > oo), but it should be remembered that is not the physical S-matrix; it still needs to be transformed 
back to the physical frame, which will remove the time-dependence (see below). ??™ qKji (i) can be evaluated explicitly 
in the WKB approximation: 



(t) = / dz 



L \ V2 

Y / U nl (z,t)(eq.K n \Q n I j\eqK n )) 

1=1 J 



where £(eqre„;i) is the classical turning point associated with the eq«„ decoupled channel, i.e., is the z value that 
solves 



Y,U nl (z,t){eqK n \Q^\eqK n ) = k 2 z . 



1=1 

To complete the solution of the scattering problem within this Matrix Diagonalization Sudden (MDS) approximation, 
as mentioned above one needs to transform back to the physical scattering matrix S" . It can be shown that this 
is accomplished by the inverse unitary transformation pi| : 

That this transformation removes the time-dependence will be demonstrated explicitly in the coordinate representation 
Sudden, below. Equations (|39|)-(|4C|) together provide the amplitude for scattering by the n th surface particle. In the 
present approach the contributions from all these particles are considered to be additive, so that the full transition 
amplitude from channel eq to e'q' is: 

N 

(e'q'|5|eq) = ^(e'q'<|5"|eqK n ). (41) 

n=l 

Note that these are generally complex-valued amplitudes, so that there may be interference between scattering from 
different surface particles. The representation of the S-matrix as a sum over particles may be questioned. However, 
since He at typical scattering energies has a wavelength between 0.5Aand 5Awhereas in a liquid the average inter- 
particle separation is at least on the order of the corresponding solid's lattice constant, i.e., several A, only at extremely 
low incidence energies this approximation is expected to break down. This however, is anyway beyond the range of 
validity of the SA. Moreover, if a representation in terms of lattice modes is required, it can be transformed to by 
representing the particle displacements in terms of creation (a^) and annihilation (a^) operators pTj: 



N 



a nt = (2m„iV^)- 1 / 2 e l e jK - R " ) 



where r1°^ denotes the equilibrium position of the n th surface atom and u>i , , and e.; denote the frequency, wave 
vector, and polarization vector of mode i, respectively J42J. 



C. Coordinate Representation Sudden (CRS) 

In practice the MDS has almost not been used; instead one can transform to the coordinate representation to obtain 
the more convenient and popular Coordinate Representation Sudden (CRS). To accomplish this it should be noted 
first that if they are not truncated, the Q nl are trivially diagonal in the coordinate representation, since they are made 
of matrix elements of the functions <5"'(R, u„,i) [^3|. Therefore, in this representation: 
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(iRu n |Q£|t'R'<) = Q nl (R,u n ,t)S(t' - t)5(R- R')5(u n - <) 



(42) 



At this point it should be clear when a single diagonalizing matrix B n , required for the MDS, can be found for all Q n 
matrices: since they are all simultaneously diagonal in the coordinate representation if they are not truncated, they 
all commute, and are therefore simultaneously diagonalizable in any other basis as well. The same reaso ning can be 
applied to the scattering matrix, which was found to be diagonal in the same basis as the Q nl [Eqs.(|37|), (|3g)]. Thus 
S n is also diagonalized by transforming to the coordinate representation: 



(iRu n |<S£|t'R'0 = S n (R, u n )S(t' - t)S(R - R')S(u n - u' n ) 



(43) 



Here ^"(R, u„) is found from the solution of the coordinate representation version of Eq.(p8|) (obtained by inserting 
the unit operator J dt' dR' du' n \t'Il'u' n ) (t'R'u' n \ before |$') and multiplying on the left by (£Ru„|): 



^2+ k D - r"(R. :.u„./i 



$lRu„(z) = 0, 



(44) 



where ^"(R, z, u„, t) was defined in Eq.(|2l|), and $tR U „ (z) — (tRu n \B n \& n ) depends parametrically on t, R and 
Eq.(p[) shows that the SA treats the coordinates t, R and u„ in an adiabatic approximation. One can now write: 



with the elastic-scattering phase shift evaluated from Eq . (|44[) in the WKB approximation as: 



77™ (R, u„ ,t) 



dz 



f(R,u„,t) 



(k 2 z - V n (R, z, u„, t)) 1/2 - k 2 z + fc z £(R, u„, t) 



Here £(R, u„, t) is the classical turning point pertaining to Eq.(44|). Next, the coordinate representation S-matrix has 
to be transformed back to the physical one, again using Eq.(fiO): 



6'q'/<JS"|6qK„) = 



J dt' dR' du' n dt dR du n (e'q / <|(B n )- 1 |tRu n ) (tRu n |<S£|i'R'u;) (t'R'u' n \B n \eqK r , 



(45) 



The integrals over the primed variables are eliminated with the help of Eq.(|43j); in addition one needs the B n matrix 
elements. On the one hand: 

QSU„, e 'q'< =<e'q'<|Q"W„) = 
(e'q'<|(B")- 1 Q2, i B n |eqK„) = 

J dt' dR' du' n dtdRdu„ (e'q'<|(S")~ 1 |tRu„)(tRu„|Q^|/;'RO(t'Ru' n |/3 ,l |eqK„) = 
J dt dR du n {(iRu n |/3"|e'q'<)}* Q" Z (R, u„, t) {(tRu n |B"|eqK n )} , 
where the last line follows from Eq.(E2j) and the unitarity of B n . On the other hand one has from Eqs.(p6|),(^9|): 



O nl , , , = 

i- fdtdR du n {/ < ,(u n )e- 1 ( e ' t / ft - ( ' , - R ) }* Q"*(R, u„, t) {/ K „,tK, t)e- l ^A-q- R ) } 



which shows that B n is nothing but the matrix of transformation between the momenta and coordinates: 

(tRu n \B n \eq Kn ) = ^-f Kn (n n ,t)e^ R e- iet / H . 
Zn 

Inserting this into Eq.(|45|) one obtains finally: 
(e'q'K' n \S n \jqK n ) = 

£^ J dtdRdu n {/:,(u„,t)e- i ^' R - £ '*/ fi '}e 2 ^( R ' u -*'{/ K „(u„,t)e i ^ R - t /' 1 )}, 



(46) 
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for the amplitude to scatter off the n th surface particle, accompanied by a momentum transfer hq — > hq' and energy 
transfer €j — > e' by the He atom, and a change TiK n — > TiK' n in the momentum of the surface particle. As in the 
MDS, the full amplitude for the scattering of the He atom is given by the sum over all surface particles, Eq. (f4l|) . The 
normalization factor a is found from the conservation of the total scattering probability: 

I de' dq' d^Ke'q'^J^bq^)! 2 = 1- (47) 

This completes the derivation of the TDSCF-SA, which is given in final form by Eqs. (£!]), © and ©. The novel 
feature in the present result is the appearance of the time dependence and the explicit accounting for the change 
undergone by the surface particles; for comparison see again the usual SA expression Eq.(||). 

The time dependence in Eq.(^6|) can be made even more explicit by expanding the wave function of the surface 
particle in terms of its energy eigenstates: 



|k„) = ^c in |j n )e" 



-%Ea t/h 



or, by multiplying on the left by (u„|, in the coordinate representation: 

/ K „(u n ;t)=^ Ci „^ n (u„)e-^ t / R , 



where <j>j (u„) = (u n |j„). Plugging this into Eq. (|46|) one obtains the matrix element: 

{K ,J e 2i V :(R,u n ,t)\ Kn) = J dUn /* ; ( Un;t ) e 2^(R.W) /Kn ( Un;t ) 



_ „/* „ j(E jL -E in )t/h 

On Jn 



^ c; k c e ' '» M jUn (R,t), (48) 



where: 



M fJn (R,i) = J ( iu n </)*,(u„)e 2 <( R ' u - i ^ in (u„). 



The Mj' n j n are the R-dependent amplitudes for a transition from energy eigenstate \j n ) to \j' n ) for the n th surface 
particle in the "£ th " configuration. The off-diagonal elements, with j' n ^ j n , correspond to inelastic transitions due 
to the collision with the He atom, as can be seen by the presence of the oscillatory term with the Ej> — Ej factor 



in Eq.(48). This demonstrates explicitly how the TDSCF-SA accounts for energy transfer from/to the liquid surface. 
The Mji^j^ have been calculated analytically (under certain simplifying assumptions) for the harmonic oscillator case 
by Schinke and Gerber ^7j, where they give rise to the Debye- Waller factor. The present formulation is more general 
as it is not restricted to phonons, and explicitly includes the energy transfer to/from the He atom. This novel feature 
should make the present result useful in a comparison with inelastic He scattering data from liquids, to be undertaken 
in a future publication. 



V. SUMMARY AND CONCLUSIONS 



This paper has introduced two novel extensions of the Sudden Approximation (SA) for atom scattering from surfaces, 
namely: (1) scattering in the presence of double collisions between the impinging atom and the surface (the Iterated 
Sudden: IS), and (2) scattering from a time-dependent quantum liquid surface (the time-dependent self consistent 
field SA: TDSCF-SA). 

The IS is a natural generalization of the SA, which takes double-collision events between the atom and the surface 
into account. An exact version was presented, suitable primarily for numerical applications, along with a more 
approximate version, which can be used for analytical work. The exact version is only slightly more computationally 
expensive than the original SA. It was shown that the IS predicts the presence of new rainbow peaks in the angular 
intensity distribution, which arise due to double collisions, and are absent in the SA. These rainbow peaks can be 
experimentally distinguished from those arising from single-collisions by their dependence on incidence energy. The 
IS seems to properly generalize the SA, which was the goal in its derivation. It is of interest to apply the IS to a 
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realistic problem, e.g. a CO adsorbate on a Pt(lll) surface, and compare the IS predictions to those of the SA and 
exact calculations. 

The TDSCF-SA generalizes the SA by taking into account the dynamics of the surface particles. The latter are 
treated as independent, with each particle moving in an average potential determined by all the others, and the 
incident atom. The resulting expression for the scattering amplitude depends on the energy transfer between the 
incident atom and the surface, and it is this inelastic contribution which constitutes the main generalization of the 
original SA approximation. The TDSCF-SA could find applications in simulations related to the novel experiments 
on inert atom scattering from molten metal liquid surfaces, such as those by Nathanson and co-workers fffi-jsif . 

Theoretical applications of the methods developed here will be undertaken in the future, but it is hoped above 
all that this work will stimulate experimentalists to further utilize inert atom scattering in the study of increasingly 
complex surface disorder. The results presented here suggest that such experiments can reveal a wealth of information 
concerning disordered surface structure, be it in the atomic-level details of isolated adsorbates, the dynamics and time- 
dependent structure of quantum liquid surfaces, or the statistics of randomly corrugated surfaces. 
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APPENDIX A: CONNECTION OF THE ITERATED SUDDEN EXPRESSIONS WITH THE 

SINGLE-COLLISION CASE 

1. Wave- Function Iterated Sudden 

It is interesting to rewrite Eq.(^) in a form which makes the connection with the single-collision case, i.e., the 
ordinary SA, more transparent. Equation (||) can be rewritten as: 

^ (kin * k ou t) = 

— ^ fdq ( dre'KKi.-qJ R-C^Jf^ f rfR ' V (R', z)e l [(q-Kout)-R'] _ ( A1) 
(2tt) 5 J J J 

The first term in the last expression is, up to a constant factor, exactly the transition amplitude from a Sudden state 
|?/f n ) to an intermediate state (q, k* ut \, as is shown next by calculating the (q, k* ut \ component of the SA wave- function: 

<q,CtlVi> = Jdx (q,Ct|r) (r|^> = 

J dq' J dr $ J q ,(z)e j(K '" +q ' ) - R e- 4(q ' fc °- ) - r (27r)"^ = 

V2^a J dr $k( z ) e *[( K --<i)- R - fe o Ut ^ 

which is indeed the term mentioned above. Having verified this, the double-collision expression obtained within the 



IS, Eq.(Al), can be interpreted as the transition amplitude from an initial state IV'fn) ( a fter a first collision) into 
an intermediate state (q, fe*„ t |, multiplied by a Born-type propagator P = J dH' V(R' , z) exp [i ((q — K out ) • R')], 
which takes the particle to the final state (k out | (after a second collision) , integrated over all intermediate diffraction 
vectors q. The "propagator" P depends on z through V(H' , z). Thus the first and second collisions are coupled not 
only by the intermediate momentum vector, but also by the (non-adiabatic) z-coordinate. As shall be shown shortly, 
this coupling is closely connected to the distinction between off-shell and on-shell amplitudes, and can be used to 
determine when an off-shell amplitude can be approximated by an on-shell one. 
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2. Phase-Shift Iterated Sudden 



Consider Eq.([lO|): the expression in parentheses is, according to Eq.(||), just (q|5 sud |Kj n ), i.e., the single collision 
amplitude for the transition from |K in ) to the intermediate state (q|, in the SA. Thus Eq.(|i"o|) can be interpreted 
as expressing the transition amplitude from an initial state K in ) to a final state |k out ), after passing through an 
intermediate state |q). The expression in parentheses in Eq.(|lO) then expresses the first collision, which is subsequently 
multiplied by a Born type propagator, expressing the second collision: the transition from (q| to k out ). But as opposed 
to the situation in the Wave- Function expression Eq.( |Al| ), the two collision events in Eq.(|lC) are decoupled in the 
z-direction. This is due to the on-shell approximation used in deriving in Eq.([[o|). Hence one would expect an on-shell 
approximation to be valid as long as the behavior of the particle in the z-dircction is primarily determined by the 
initial z-component of the particle's wave-vector. For in this case, the influence of the potential on the propagation 
in the z-direction is relatively negligible. This is in the spirit of the high-energy SA. 
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FIG. 1. On the left, the first collision strikes the adsorbate, the second strikes the surface. The reverse sequence occurs on 
the right. 



